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I. INTRODUCTION 

The ternary intermetallic system Al-Pd-Mn has been 
of great interest in the last years, because it forms a high 
number of complex metallic alloy compounds (CMAs). 
In this paper we focus on the S phases, which are approx- 
imants of a decagonal quasicrystal with a lattice constant 
of 1.6 nm in the periodic direction. Under plastic defor- 
mation, these phases show a novel type of dislocations, 
so-called metadislocations, which were first described by 
Klein et al^ 

Ab initio studies of these metadislocations, even w ith 
fast codes using density functional theory like VASP,!^ 
are currently unfeasible. Their spatial extent is about 
200 A and they involve more than 10 000 atoms - im- 
possible to simulate even with state of the art ab initio 
programs. 

With classical molecular dynamics (MD) it is easily 
possible to simulate structures with millions of atoms in 
reasonable time. The treatment of atoms as point masses 
interacting with an effective potential allows for micro- 
scopic insight into many processes on the atomic scale. 
The ability to control almost any aspect of the simula- 
tion can be used for optimizing the structure, determin- 
ing physical properties or explaining physical phenomena 
in detail. 

However, obtaining an effective potential for classical 
molecular dynamics is not straightforward. In order to 
extract reliable results, a potential has to be adjusted 
to the specific physical conditions considered. These can 
be for example high pressures, strain, surfaces or phase 
boundaries. A common way is to fit a potential such that 
it reproduces experimental dat a like lattice constants, co- 
hesive and surface energied^ or simply combining pure 
element potentials into an alloy potential. 

For ternary systems, like Al-Pd-Mn, establishing a 
potential with these approaches is very challenging. 
The small number of available experimental data is 
not enough to fit reliable effective potentials. Hence, 
to obtain a potential that can be used for structure 
analysis and optimization, we apply the force-matching 



method^ using the potfit package.^ * In the force- 
matching method, results from ab initio simulations are 
used as reference data to adjust the parameters of a po- 
tential. This not only dramatically increases the amount 
of information available for fitting (the total number of 
datapoints can easily reach several thousands). Also, if 
the reference data is found to be insufficient, more per- 
tinent reference data can be generated at relatively low 
cost. This makes it possible to create realistic potentials 
for binary or ternary systems. In our case we used forces 
on individual atoms, the cohesive energy and stresses on 
the unit cells to fit a reliable potential. 

In Sec. In] we describe the interaction model used in this 
research. The fitting procedure using the force-matching 
method is presented in Sec. |111[ the reference data used 



is given in Sec. IV The results will be discussed in detail 
in Sec. |Vl 



II. EAM POTENTIALS 

A common way to describe atomic interactions in met- 
als is the embedded atom method (EAM).'2llt implicitly 
includes many-body interactions by a term which de- 
pends on the environment of every atom. The potential 
energy of a system described with the EAM method can 
be written as 



(1) 



with ni^'^pj{rij). (2) 
The first term in ([l]) represents the pair interactions be- 



tween atoms i and j at a distance r.^ 



The 



function Fi{ni) is the embedding energy of atom i in the 
host density n^. This density ^ is calculated as the 
sum over contributions from the neighboring atoms, with 
Pj being the transfer function of atom j. It does not rep- 
resent an actual physical density; rii is a purely empirical 
quantity. 
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For the pair and transfer part, we have tested three 
different combinations of analytic functions as model po- 
tentials. Potential I has oscillations in the pair potential 
but not in the transfer function. In contrast, potential II 
has oscillations only in the transfer function. Finally a 
third potential has oscillations in both functions. 

For the simple pair potential without oscillations we 
chose a Morse potential. It has a single minimum and is 
used in model II only: 
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^' is a cutoff function, where the free parameters Vc and 
h describe the cutoff radius and the smoothing of the po- 
tential. The remaining parameters are De, a, and r^; 
is the depth of the potential minimum, r^, the equilibrium 
distance and a the width of the potential minimum. The 
pair potential function with oscillations is adopted from 
Mihalkovic et al.W 



<i>(r) = * 



cosl/er 4- ip) 



(4) 



This "empirical oscillating pair potential" (EOFF) has 
been used in various works on complex metallic alloys and 
quasicrystalSj'^i^^ as it provides great flexibility. The 
first term of Q with the parameters C\ and r^i controls 
the short-range repulsion. The second term is responsi- 
ble for the damping (C^tTi-i) of the oscillations with the 
frequency k. 

The cutoff function "^(x) is defined by 



X 



(5) 



for X < and "^(x) = for x > 0. This function guaran- 
tees that the potential functions as well as their deriva- 
tives up to the second order approach zero smoothly at 
the cutoff distance r^. 

Two different analytic forms were used as transfer func- 
tions; one allows for oscillations, the other one does not. 
The latter one is a simple exponen tial de cay frequently 
used in established EAM potentialsPI^ 



p(r) = Q;exp(— /?r). 



(6) 



where a is the amplitude and /? is the decay constant. 
This function is used in model I. For models II and III, 
we used an oscillating transfer function, which is taken 
from Rcf. [H 



p{r) = * 



r — r. 



1 + ai cos(ar) -|- a-i sin(ar) 
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The four free parameters are ai,a2,a and /3, where oi 
and 02 determine the amplitude of the oscillations, a is 
the wave vector and /3 controls the decay. 

The embedding function F{n) was adopted from 
Ref. It is based on the general equation of state from 



Rose ei ai^^ The original form is given as 
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The parameters in this function are Fq, Fi,p,q and rig. 
p and q are real values and rig is the equilibrium density. 
In this paper we use this function in the limit p ^ q and 

chose Hp = 1: 



F{n) = Fo [1 - glogn] n'^ + Fm, 



(9) 



because the original form is numerically unstable with 
our optimization algorithms. 

The number of free parameters of our three potential 
models is comparatively large. The non-oscillating (oscil- 
lating) pair potential has 3 (6) parameters, and the non- 
oscillating (oscillating) transfer function requires 2 (4) 
values. All models share the embedding function with 
3 free parameters. Every pair and transfer function has 
one additional parameter h for the cutoff function 5*. The 
cutoff radius is kept fixed at 7 A. In a ternary system 
like Al-Fd-Mn with 12 potential functions, this adds up 
to a total number of 60, 48 and 66 parameters for the 
models I, II and III, respectively. 



III. FITTING PROCEDURE 

All force-matching was performed with the potfit 
package of Brommcr and G abler,''''®' which has previously 
been used to optimize tabulated pair and EAM poten- 
tials. For this work, its capabilities were extended to 
analytic potential models. 

All free parameters of the analytic functions were fit- 
ted to an ab initio reference database containing relaxed 
(T = 0) structures, snapshots from ab initio MD simu- 
lations at higher temperatures and a few strained sam- 
ples (see Tables [l| and [ll| . All ab initio calculations were 
performed with the Vienna Ab initio Simulation Pack- 
age (VASP)2EI using the generalized gradient approxima- 
tion (GGA) and the Projector Augmented Wave (PAW) 
method.^® 

Two different optimization algorithms were used to fit 
the potentials. They both minimize the sum of squares 
defined by 



(10) 



where AE, AF and AS" are the energy, force and stress 
residuals. These deviations are calculated as the differ- 
ence of the ab initio and the EAM value, e.g. 



AE = E 



EAM 



E, 



ab initio- 



lje and ws are global weights for the energies and 
stresses, loe — 22 500 was chosen to obtain poten- 
tials that yield very precise energies, but also reason- 
able forces. For configurations with about 150 atoms. 
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this effectively weiglis tfie energies witli a factor of ap- 
proximately 50. The stress weight ujs was set to 750, so 
that the total weight of the six stress tensor components 
per configuration is approximately equal to ten times the 
weight of all forces in one configuration. 

The first optimization algorithm used is simulated 
annealing It is based on the Metropolis criterion, where 
a decrease in the target function Z is always accepted 
and an increase only with a probability P ~ e^^^^'^ . 
This allows the algorithm to escape local minima. The 
artificial temperature T is steadily decreased during the 
optimization. To ensure that the fit converged to the 
global minimum, the optimization was restarted with a 
high temperature several times. Subsequently a conju- 
gate gradient based method^ was applied to converge 
to the final optimum. During the fitting procedure, all 
parameters were confined to a predefined range by use of 
numerical punishments. 



IV. REFERENCE DATA 

The structures used as reference data are shown in 
Tables |T] and [n] There are 119 configurations with a total 
of 16 103 atoms. The number of reference datapoints is 
49 340. They consist of 48 309 forces, 119 energies and 
714 stresses. 



Number of atoms 



Ai? (cV/atom) 



AI-Mn structures 


Al-Pd structures 


AlioMn3./iP26 


AlPd.cPS 


AliiMn4.aP15 


AlaiPdg.t/lie 


Ali2Mn.c/26 


AlaPda./iPS 


Al6Mn.oC28 




KlMn.tPA 





TABLE I. Binary structures (T = 0) used to fit the potentials, 
with their corresponding Pearson symbol. 

In addition to the binary and ternary structures, one 
reference configuration for each of the pure elements was 
also included. These were, in detail, Al.cF4, Pd.ci^4 and 
Mn.c/58. This was done to get reliable reference points 
for the calculation of the enthalpy of formation. 

All atomic configurations from binary systems (Ta- 
blejl]) were taken from the alloy database of Widom et 
aZ.l^and have been fully relaxed with ah initio meth- 
ods. They were chosen to provide more data for the 
Pd-Pd and Mn-Mn interactions. Magnetism was not 
included in our ah initio calculation; it was shown that 
the manganese atoms in the S-phases are nonmagnetic.^'^ 
Because the structures we want to investigate are on the 
aluminum-rich side of the phase diagram, there is only 
little data for the Mn-Pd interaction. 

The reference configurations for the S-phases are from 
different sources. The structures in Tablelllldenoted with 
superscript a were taken from the alloy database.'^^lThey 
were generated with the canonical cell tiling ,1^3 which cre- 
ates hypothetical models by decorating a tiling with clus- 



T = AlgaPdagMmd 
Al92Pd28MnJ 
Alii2Pd36Mn| 
Alii4Pd34Mne^ 
Alii2Pd34Mnr 
Aliio+,Pd32Mn£] 

Ali24Pd8Mn2ffl^ 
Ali47Pd43Mnis* 

Al294Pd88Mni£'' 

r>0 Al92Pd28MnflT 

Al92Pd28Mni([2](Mro K) 



-0.512 
-0.485 
-0.526 
-0.503 
-0.512 
see Sec 



VA 



-0.297 
-0.485 
-0.491 



^ Structure generated from canonical cell tilinj^ 
^ From structure optim izat ion 
" T-Al-Pd-Mn, see Sec. |IV] 

From several MD runs at 600, 1100 and 1800 K with small 

strains 

TABLE II. Ternary structures used to fit the potentials and 
their ah initio formation enthalpy A_H". 



ters. To compensate for the low amount of manganese in 
these samples and the hence resulting lack of data, five 
of the aluminum atoms were replaced by manganese in 
some configuratio ns. Ah initio molecular dynamics simu- 
lations with VASFf2ini were run with these samples at 600, 
1100 and 1800 K to obtain different local atomic configu- 
rations. These calculations were done in the generalized 
gradient approximation (GGA) with PAW potentials. ^^1 

At the same time, ah initio structure optimization was 
carried out for two of the 5-phases. Particularly this 
were the S-phase with the smallest unit cell, which con- 
tains about 152 atoms and is called ^ and the next bigger 
one, containing about 304 atoms, which is called All 
structures generated in the course of this optimization 
are denoted in Table [TT| by superscript b. 

To judge the stability of these structures, their energy 
is compared to a mixture of competing phases, the con- 
vex hull. This hull, defined over a ternary phase diagram, 
contains the cohesive energies of all stable compounds as 
vertices. If the energy of a structure is above this convex 
hull, it could decompose into the neighboring structures 
and thus lower its energy. If the energy of a new structure 
is below the convex hull, it is considered to be thermody- 
namically stable. The structures which define the convex 
hull for the S-phases, are T-AlPdMn, Ali2Mn, Al2iPd8 
and Al3Pd2. They have also been included in the ref- 
erence database. A detailed description of these phases 
and the convex hull is given in Ref. 



V. RESULTS 

We determined parameters for all three potential 
models from the reference data described above. The 
root mean square (RMS) errors for forces, energies and 
stresses after the optimization are in the same order of 



magnitude for all models (see Table III I . While model III 
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has the smallest errors for forces and energies, model I 
has the biggest errors for all three quantities. Model II 
has the smallest stress deviations. While the force error 
for model I is about 20% larger than the one for model 
III, the energy error is significantly larger with about 50% 
diff'erence. 



RMS errors for 


Model I 


Model II 


Model III 


forces 


265.63 


221.40 


220.07 


energies 


19.36 


14.49 


12.53 


stresses 


99.99 


76.83 


98.30 



TABLE III. Root mean square errors after the optimization 
for forces (in meV/A), energies (in meV/atom) and stresses 
(in kPa). This data is calculated with the reference configu- 
rations used for fitting the potentials. 

A graphical representation of these errors can be seen 
in Fig. [T] The scatterplots in the upper row display the 
energies of the reference data. Forces are shown in the 
lower row. The range of the force plots is due to the many 
high temperature MD simulations that are included in 
the reference data. The forces therein can become very 
large because of the short interatomic distances that may 
occur at these temperatures. 



Model I Model II Model III 




-6 -3 3 6 -6 -3 3 6 -6 -3 3 6 



FIG. 1. Scatterplot for energies and forces with the EAM 
values on the vertical axis and the ab initio reference data on 
the horizontal axis. The insets are magnified by a factor 4.5. 

These errors cannot solely be used to judge the quality 
and transferability of the potentials. For that purpose 
another set of ab initio data has been extracted from the 
structure optimization. It has not been included in the 
reference data and can be used to determine the trans- 
ferability of the different potentials. The same errors as 
before have been calculated and can be seen in Table 
|IV[ As with the reference data, model III has the lowest 
force and energy errors. The relative error of the en- 
ergy is about 0.2%, for stresses about 5% and 550% for 
forces. This is due to the fact that all configurations in 
this test data are ground state structures and therefore 



only contain very small forces. 



RMS errors for 


Model I 


Model II 


Model III 


forces 


141.90 


131.90 


130.46 


energies 


10.42 


10.47 


10.28 


stresses 


32.39 


23.76 


36.89 



TABLE IV. Root mean square errors for forces (in meV/A), 
energies (in meV/atom) and stresses (in kPa). This data is 
calculated with test data, containing only structures that were 
not included in the optimization process. 

The errors for the test data in Table IIVI are smaller 
than those of the reference configurations in Table |III[ 
because there are only ground states included and no 
high temperature MD runs. 

Based upon these simple energy and force considera- 
tions, all the potential models appear to be of similar 
quality. Model III, however, should be slightly superior 
to the other two potentials. Further tests are necessary to 
determine the performance of the potentials in different 
situations. They will be presented in Subsection |V B| 

A. Structure Refinement 

In Ref. "M", the structure of the S-phases of Al-Pd-Mn 
has been optimized by energy minimization in ab initio 
and molecular dynamics simulations. We use several of 
the structures tested there to judge the quality of the 
optimized potentials. The S-phases consist of columns 
of pseudo-Mackay icosahedral clusters (PMIsVP^ a slight 
deviation of the famous Mackay icosahedron.l^Sl 

Every PMI cluster consists of a single atom at the cen- 
ter with a first shell of an experimentally poorly deter- 
mined number of aluminum atoms. The second shell is an 
icosahedron of 12 transition metal atoms and the outer 
shell an icosidodecahedron of 30 aluminum atoms, see 
Fig. [2] Almost all atoms of the S-phases belong to these 
clusters. It is difficult to measure the exact number of 
atoms in the first shell because aluminum atoms are hard 
to observe in diffraction experiments. 

For the lowest quasicrystal approximant, the ^-phase, 
there are four PMI clusters in one unit cell. Several dif- 
ferent occupancies of aluminum atoms in the first shell 
were tested in Ref. 24, Each configuration was denoted 
by a single number, giving the average number of alu- 
minum atoms per cluster. Structures from eight up to 
eleven atoms per PMI were generated and tested. 

The results with the different potential models can be 
seen in Table |Vj All structures were completely relaxed 
with ab initio methods, the corresponding ab initio en- 
ergy is given in column 2. The energies of these config- 
urations with the generated EAM potentials have been 
calculated after subsequent relaxation with the respective 
potentials. This relaxation causes small displacements of 
the atoms from their ab initio determined positions. For 
models I these average displacements are O.IO A/atom, 
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FIG. 2. (color online) Detailed structure of the pseudo- 
Mackay icosahedral cluster. A few atoms of the second and 
third shell are omitted to see the central atom and the alu- 
minum atoms in the first shell. The icosahedron and icosi- 
dodecahedron are indicated by planes and bars, respectively. 
The central atom is manganese. Aluminum atoms are de- 
picted in yellow (light gray), and palladium atoms in magenta 
(dark gray). 



0.08 A/atom for model II and 0.11 A/atom for model III. 
This clearly shows that all potential models can stabiUze 
the ground states of all structures that were generated. 



Number of atoms 
per PMI 


Eab initio 

(eV/atom) 


AE (meV/atom) 
Model I Model II Model III 


8 


-4.753 


-13 


-12 


-20 


8.25 


-4.755 


-6 


-7 


-13 


8.5 


-4.756 


-1 


-3 


-5 


8.75 


-4.757 


+3 


+2 


+2 


9 


-4.755 


4-4 


+4 


+3 


9.25 


-4.747 


+ 1 


+ 1 


+1 


9.5 


-4.741 


+ 1 


+0 


+2 


9.75 


-4.731 


-6 


-5 


-2 


10 


-4.731 





+2 


+4 


10.25 


-4.714 


-12 


-12 


-5 


10.5 


-4.704 


-15 


-17 


-7 


10.75 


-4.692 


-19 


-21 


-13 


11 


-4.683 


-22 


-24 


-17 



TABLE V. Cohesive energies (in eV/atom) of different opti- 
mized configurations for the ^-phase. The energy differences 
AE between the ab initio calculations and the respective 
model are given in meV/atom. 

All models are having difhculties with the energies of 
structures that contain less than 9 or more than 10 atoms 
in the inner shells of the PMI clusters. This may be 
an indication for the mechanical instability found dur- 
ing the structure optimization.'^ The energy of these 
structures is highly unfavorable; at elevated tempera- 
tures some atoms drifted from the outer shell to the inner 



shell or vice versa to achieve an inner shell with 9 or 10 
aluminum atoms. 

All energy differences between the ab initio and EAM 
calculations are smaller than 10 mcV/atom for config- 
urations ranging from 8.5 to 10 atoms per PMI cluster. 
This energy is considered a critical threshold for the accu- 
racy of the potentials. Regarding the energy differences 
between the different structures, which are on the order 
of 1 meV/atom, all potentials can evidently distinguish 
between these different configurations. 

The structure optimization in Ref . [211 yielded four al- 
most stable structures, which are different from the ones 
shown in Table |V) There, not only the atoms in the in- 
ner shell are varied, but also atoms not belonging to the 
PMI clusters. These alterations were not done in a sys- 
tematic manner, the structures will be listed in tabular 
form. The amount of atoms for ^- and ^'-phases is the 
same, only the arrangement of the PMI cluster columns 
is different. These structures were tested with the three 
different potentials. The results can be seen in Table VI 
The two upper structures in this table are ^-phase 
lower two structures are 



the 



Ea 



AE (meV/atom) 



composition 


(eV/atom) 


Model I 


Model II 


Model III 


C-228-64-12 


-4.702 


-5 


+4 





C-224-68-12 


-4.748 


+1 


+7 


+4 


^'-228-64-12 


-4.703 


-5 


+3 


+ 1 


^'-224-68-12 


-4.748 


+1 


+5 


+5 



TABLE VI. Cohesive energies (in eV/atom) of the four almost 
stable phases after relaxation. The energy differences A_E are 
given in meV/atom. The composition is given in numbers 
of aluminum, palladium and manganese atoms, in this order. 
All configurations have 9 aluminum atoms in the inner shell 
of the PMI clusters. 

After the relaxation with the effective potentials, all 
models show a very good agreement with the ab initio 
calculated energies. The mean displacements after the 
relaxation are again in the same order of magnitude as 
before, 0.11 A/atom for model I, 0.08 A/atom for model 
II and 0.15 A/atom for model III. 

Based on these pure energy comparisons, all three po- 
tential models seem to be of equal quality, with slight 
advantages for model III. 



B. 



Tests 



A force-matched potential is only useful, if it can repro- 
duce key quantities that were not directly included in the 
reference data. Here, we subjected the three potentials 
to a series of tests. The first test is whether the poten- 
tial can stabilize the ^-phase even at elevated tempera- 
tures. As there was a large number of high temperature 
ab initio MD simulations included in the optimization, 
the potentials should be able to preserve the structure 
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of the ^-phase under these conditions. We carried out 
an ab initio MD simulation at 1200 K for 50 ps)^ where 
the phase is still mechanically stable. In a time-averaged 
picture of the density, the atoms in the two outer shells 
of the PMI clusters did not move, but the atoms in the 
first shell showed some rotational degree of freedom. 

All three models were able to stabilize the structure 
at this temperature. While models II and III give the 
same results as the ab initio calculation (cf. Ref. 
model I shows additional degrees of freedom. In the time- 
averaged picture the atoms forming the outer shell of the 
PMIs are not as steady as in the ab initio simulation. 
Also the atoms, which do not belong the these clusters, 
exhibit a density distribution that is twice as large as 
expected. This means that model I may have difficulties 
stabilizing the structure at even higher temperatures or 
against fluctuations in the local atomic arrangement. 

For molecular dynamics simulations, the stabilization 
of different phases can be a problem. We checked some 
well known phases for all three potential models with re- 
spect to cohesive energy and phase stability. The results 
can be seen in Table |VII[ All three potentials can sta- 
bilize the different phases. The deviation of the atomic 
positions after relaxation compared to the ab initio ref- 
erence values is very small. The energies are reproduced 
with errors of under 200 meV/atom. 

Another important test is the calculation of formation 
enthalpies AH with the potentials. AH is defined as the 
energy difference of a structure to the tie plane of the 
pure element energies. This has been calculated for all 
configurations in Tables [V] and | VI[ The reference energies 
are given in Table [VlTTj For the structures with different 
amounts of aluminum atoms in the inner shell of the PMI 
clusters, the results can be seen in Figure [3] The devi- 
ations from the ab initio enthalpies are very similar to 
those from Table fVl For less than 8.5 and more than 10 
atoms in the inner shell of the PMI clusters the enthalpies 
differ more than 10 meV/atom. 

The enthalpies for the four almost stable structures are 



shown in Table IX All three models give very accurate 
enthalpies with deviations all smaller than 10 meV/atom. 

During the structure optimization a very long ab ini- 
tio MD run with 50 000 steps at 1200 K was performed. 
Snapshots were taken from this simulation at different 
timesteps and quenched very rapidly. This has also been 
done with the EAM potentials. The results show a very 
good agreement for different snapshots. The structures 
only differ very slightly in atomic positions. While there 
is a steady offset of about 100 meV/atom in the energy 
for higher temperatures, the overall trend can clearly be 
followed. For lower temperatures and T = the energies 
were in the same order as for the structures in Table [V] 
There were no major differences for all three potential 
models. 

To determine if a structure is thermodynamically sta- 
ble, the energy difference of this structure to the convex 
hull is calculated. If this difference is negative, the struc- 
ture is stable, otherwise it it unstable. For more details 
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FIG. 3. Comparison of the ab initio formation enthalpy AH 
(in eV/atom) with the three potential models. The lines be- 
tween the ab initio datapoints are added as a guide to the 
eye. 



on the convex hull see Ref. [24l This energy difference 
has been calculated for all structures in Table |V] and is 
shown in Fig. |4] 



45 
40 
35 
30 
25 
20 
15 
10 
5 







ab initio 














Model I 


--A-- 








Va 




Model II 


--G-- 








" \ "A 




Model III 


■■-B-- 


























~&---A, 












''G-- 






-_A-- 


A- 


-A 










-o- 


--0 - 


8 8.5 


9 


9.5 


10 


10.5 




11 



average number of atoms per PMI cluster 

FIG. 4. Difference of the cohesive energy to the convex hull 
for different amounts of aluminum atoms in the inner shell of 
the PMI clusters. The datapoints for model I are shifted by 
-|-10 meV/atom and by -|-5 meV/atom for model II. 



For the sake of clarity, the datapoints of model I and 
II are shifted by -1-10 and -1-5 meV/atom. While these 
models show a clear decrease of the energy difference 
with increasing number of atoms inside the PMI clus- 
ter, model III has minima for 9 and 10 atoms, like the ab 
initio reference calculation. As this is the main criterion 
for performing a structure optimization, models I and II 
cannot be used for this purpose. Only model III is able 
to reproduce the shape of the ab initio calculation. 

The melting point for the ^-phase has been determined 
with all three potential models. In MD simulations the 
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Eab initio 




Model 1 






Model II 






Model III 




System 


E (eV/atom) 


-Eeam 


AE 


As 


-Eeam 


AE 


Ax 


-Beam 


A-E 


Ax 


AlPd.cP2 (B2) 


-5.330 


-5.430 


-0.100 





-5.503 


-0.173 





-5.445 


-0.115 





A1Pdci rF^6 CnOt) 


-5.236 


-5.426 


-0.190 





-5.424 


-0.188 





-5.430 


-0.194 





AlaPd.i/S (DO22) 


-4.421 


-4.546 


-0.125 





-4.540 


-0.119 





-4.560 


-0.139 





Al3Pd.cP4 (LI2) 


-4.609 


-4.647 


-0.038 


0.13 


-4.651 


-0.042 


0.17 


-4.650 


-0.041 


0.17 


AhUn.tlS (DO22) 


-5.132 


-5.053 


0.079 


0.04 


-5.129 


0.003 


0.01 


-5.175 


-0.043 





Al3Mn.cP4 (LI2) 


-5.032 


-5.187 


-0.155 





-5.180 


-0.148 





-5.197 


-0.165 






TABLE VII. Cohesive energies for different phases in the Al-Pd-Mn system. All energies and energy differences are given in 
eV/atom. The mean square displacements (Ax) after relaxation are given in A/atom. A displacement of means the value is 
smaller than 10~* A/atom. 



Eab initio A-E (meV/atom) 

(eV/atom) Model I Model II Model III 



A1.CE4 


-3.688 


-5 


-4 


-3 


Pd.cE4 


-5.199 











Mn.c/58 


-8.964 


+68 


-7 


+3 



TABLE VIII. Ab initio energies (in eV/atom) and the differ- 
ences for the effective potentials (in meV/atom) for the pure 
elements. These energies were used for the calculation of the 
formation enthalpies AH. 





AHah initio 


A-ffsAM 


-AH, 


lb initio 


composition 


(eV/atom) 


Model I Model II 


Model III 


C-228-64-12 


-0.488 


-6 


+1 


-1 


C-224-68-12 


-0.513 





-H4 


+3 


^'-228-64-12 


-0.488 


-6 








^'-224-68-12 


-0.514 





+2 


+3 



TABLE IX. Ab initio formation enthalpies AH (in eV/atom) 
of the four almost stable phases and the differences for the ef- 
fective potentials (in meV/atom) after relaxation. The com- 
position is given in numbers of aluminum, palladium and 
manganese atoms, in this order. All configurations have 9 
aluminum atoms inside the PMI clusters. 



volume per atom has been calculated while the sample 
was heated from 950 K to 1400 K. At the melting point, 
there is a jump in atomic volume, which corresponds to 
the melting transition. For model I this was found at 
1130 K, for model II at 1370 K and for model III at 
1300 K. With this method, the melting point is generally 
overestimated, due to the high heating rates. For the 
simulations we chose a heating rate of 5 x 10^^ K per 
timestep, this equals 5 x 10^ K/s. If one compares these 
temperatures with the experimental value of 1118 K, the 
value for the potential model I seem to be too low, model 
II and III are in the expected temperature range. 

Another test we performed is the calculation of the 
elastic constants. All S-phases have an orthorhombic 
unit cell. The corresponding nine elastic constants were 
determined by examining the cohesive energy during ho- 
mogeneous deformations of the sample. 

The results with all three models (see Table |x| show 





ab imtio 


Model I 


Model II 


Model III 


Cii 


175.79 


255.25 


244.66 


200.98 


C22 


192.75 


269.79 


246.74 


193.61 


C33 


227.46 


243.53 


246.64 


160.57 


C12 


58.76 


158.83 


145.57 


102.76 


Cl3 


67.85 


146.75 


146.78 


92.95 


C23 


56.34 


151.19 


146.51 


107.04 


C44 


72.54 


42.57 


42.42 


42.77 


C55 


67.77 


41.46 


47.19 


46.66 


Cee 


71.25 


48.51 


48.21 


43.76 



TABLE X. Elastic constants of ^-Al-Pd-Mn in GPa. 



only very little agreement with the ab-initio values. Only 
model III can reproduce Cn and C22- AH other elastic 
constants differ by up to a factor of 3. The potentials 
are apparently not able to reproduce the shear stress. 
However, this behavior is to be expected, if one takes 
into account that these potentials were generated for en- 
ergy minimization purposes. For other applications, like 
calculating mechanical properties, an extended database, 
containing enough data on shears, should be used. The 
only samples used for these potentials, that included de- 
formations, were high temperature ab-initio MD snap- 
shots. These were strained along either of the cartesian 
axes, which are perpendicular to the periodic stacking 
axis of the quasicrystal. The corresponding elastic con- 
stants are Cn and C22, which are the only ones correctly 
reproduced by model III. 

This clearly shows that force matched potentials are 
limited in their applications. They give very accurate 
results regarding the energy and forces because they are 
tuned to these quantities. For other physical properties, 
like elastic constants, the potentials are less accurate. 



VI. SUMMARY 

The Al-Pd-Mn potentials presented are very well 
suited to model the energetics of the S-phases. They 
were obtained with the force-matching method, which is 
fitting the parameters to a large database of ab initio 
determined reference data. All three analytic potential 
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models tested were able to reproduce the ab initio values 
of the energies with very high accuracy. The error sum of 
the fitting process for all three potentials is very similar, 
yet they show very different properties when used in MD 
simulations. 

The differences of the models become visible when cal- 
culating energy differences like formation enthalpies or 
the convex hull. There, model III shows the smallest de- 
viations and can reproduce the ab initio values with very 
high accuracy. The models I and II also give very good 
energies differences but cannot be used to predict the 
stability of a structure with the calculation of the convex 
hull. This indicates that oscillations on two length scales, 
like in model III, are necessary. However, the reasons 
for this are unclear. For further structure determination 



and analysis of the metadislocations in the ^-phases, the 
model potential III will be used. 
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EOPP pair function 



pair 


Ci 




C2 




k 




h 


Al-Al 




7.6769 


—0 0333 


1.0012 


3.7658 


3.8484 


1.3897 


AI-Mn 


338.7250 


7.5484 


-0.4212 


1.9271 


2.7530 


0.0033 


0.5000 


Al-Pd 


981.8107 


9.1908 


—89.9193 


4.7322 


0.2491 


1.3235 


0.6211 


Mn-Mn 


3.8460 


19.9995 


—44.5953 


4.1469 


1.2084 


1.0115 


1.5938 


Mn-Pd 


12.8931 


3.4348 


-90.3824 


4.4851 


1.6212 


0.0005 


0.5007 


Pd-Pd 


6625.3081 


9.5962 


99.8792 


6.1164 


3.8088 


2.5086 


0.5235 








transfer function 










element 


CLl 




a 


B 


h 






Al 


0.1317 


0.0399 


2.7507 


2.3142 


1.9995 






Mn 


-1.5432 


1.0321 


1.6018 


2.4154 


1.9996 






Pd 


0.4962 


0.7317 


2.9972 


3.4308 


0.5001 












embedding function 










element 




Fi 












Al 


-2.9403 


0.5639 


-1.3026 










Mn 


-1.5862 


1.3917 


-5.3935 










Pd 


-4.0016 


0.9432 


-5.7749 











cutoff radius 
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TABLE XI. Parameters for the model III EAM potential with r in units of A and V(t) in eV. 
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FIG. 5. (color online) Plots of the 12 functions of the EAM potential (model III) for Al-Pd-Mn. 



